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Abstract 

To obtain a light mode in two-dimensional staggered fermions, we introduce four new 
local operators keeping the rotational invariance for a staggered Dirac operator. To split 
masses of tastes, three cases are considered. The mass matrix and the propagator for free 
theories are analyzed. We find that one of three cases is a good candidate for taking a 
single mode by the mass splitting. In the case, it is possible that a heavy mode obtains 
infinite mass on even sites or odd sites. 



* E-mail: 
t E-mail: 
E-mail: 



hatake@muse.sc.niigata-u.ac.jp 
hide@muse . so . niigat a-u .ac.jp 
so@muse.sc.niigata-u.ac.jp 



1 Introduction 



Staggered fermions are formulated in which species doublers of a Dirac field are interpreted 
as physical degrees of freedom, tastes, on lattice Although the fermion determinant has 

many advantages in the cost of its numerical calculation [31 13 5 it remains for a 4-fold degeneracy 
problem of tastes in four dimensions to be unsolved. A fourth-root trick of the determinant in a 
staggered Dirac operator is an approach to unfold the degeneracy and studies on its theoretical 
basis are developed 0101 EI- However, we have no local expression of one taste Dirac fermion 
after the fourth-root trick. 

Avoiding the trick, there are pioneering works for solving the degeneracy tried by improved 
staggered fermion approaches [Hlin]- The improved actions generally include more operators 
than the original staggered one and are difficult to treat them ^Uj. For the control of their 
operators, we make use of staggered fermions on a D-dimensional lattice space based on an 
SO {2D) Clifford algebra ^T]. The formulation by the S0{2D) Clifford algebra is powerful in the 
control of fermion operators and we can describe any improved fermion action on a hypercubic 
lattice. In addition, a discrete rotational symmetry (cubic symmetry) can be represented by 
the algebra. 

In this article, we analyze the mass splitting of degenerate tastes by adding four operators 
to the original staggered action in two dimensions. Only these four operators keep the discrete 
rotational symmetry in any dimension ^Tj^. The total mass matrix analysis is insufficient 
because the matrix does not commute with the kinetic term. Therefore, we also analyze the 
propagator and the pole. It is found that only one combination in these operators is a good 
candidate after these analyses. 

This article is organized as follows. In section 2, staggered fermions are formulated by the 
S0{2D) Clifford algebra. Four operators are introduced to obtain taste-splitting masses. These 
operators keep the discrete rotational invariance. We analyze the mass matrix and the mass 
pole of the improved free staggered Dirac operator in sections 3 and 4. Section 5 is devoted to 
further analyses of the massless limit and infinite mass of a heavy mode. We summarize and 
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Figure 1: A two-dimensional lattice unit and the weight of a spinor representation in 5*0(4). 
discuss about our approach in section 6. 

2 Formulation of Staggered Fermions and Rotational Sym- 
metry 

The formulation of staggered fermions on the D-dimensional lattice space has been presented 
based on the SO {2D) Clifford algebra ^T]. The basic idea is that the dimension of the total 
representation space including spinor and taste spaces, 2^ is the same as that of an SO {2D) 
spinor representation. 2^ is also the same as the number of sites in a D-dimensional hypercube. 
To avoid the double counting of sites, the lattice coordinate is noted by 

n^ = 2N^ + c^ + r^, (1) 

where is the global coordinate of the hypercube. In this case, a fundamental unit is 2a, 
where a is a lattice constant, and is set to unity, = 1/2 for any fi means the coordinate of 
a center in the D-dimensional hypercube and does the relative coordinate of a site to the 
center. The relative coordinate is the same as a weight of the spinor representation in SO {2D). 

Although our formulation can be generalized, we consider a theory in a two-dimensional 
lattice, for simplicity. Relative coordinates of four sites around a plaquette are written by 

(n, rs) = (-1/2, -1/2), (-1/2, 1/2), (1/2, -1/2), (1/2, 1/2), (2) 
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as shown in Figure Actually, our staggered fermion is defined on sites (j2I) as 



^(n) = ^UriN) 
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(3) 



\ ^(1/2,1/2) / 

It is noted that \l'i and \E'4 are put on even sites and \E'2 and \E'3 are put on odd sites. 

An 5*0(4) Clifford algebra plays a crucial role in two-dimensional cubic lattice formula- 
tions pi] . The original staggered fermion action |2] can be written as 



N,N' ,r,r',fi,T 



(4) 



where r is a two-dimensional vector with its components of ±1/2 and DJ^ for fi = 1,2, is a 
generalized difference operator defined by 

1 



" CT=0,1 



with 



(^f, = 0, 



(5) 



(6) 



a is a two-dimensional vector dual to r and (V~) implies a forward (backward) difference 
operator along the /i-direction, respectively. In Eq. ©, a link variable U2n+3,h is introduced 
for gauge covariance. This formulation of staggered fermions has usual gauge interactions and 
definite 0(a) higher order terms. The matrix F^^^ in our action is composed of the 5*0(4) 
Clifford algebra 7^ and 7^ 
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72, 



(7) 



= 7l5 -^r2,(_l/2,l/2) 



/ -i 0\ 

z 

i 

V -i J 



= 72, 



F^,.^ ^ (z7i7i)'/'+^Hn272)^/^+^^7M. 



l/2+r2. 



(9) 



Here we denote the fundamental algebra, or the SO {4) Clifford algebra as 



{7/., 7u} = {7m' %} = 25^'^, 



{7/x, %} = 0, 



(10) 



and 



{r^,,^,r5} = o, (12) 

where Fs = 71727172 = diag{l, —1, —1, 1). Prom the algebra, we find that a massless staggered 



fermion has an even-odd symmetry 



erg 



(13) 



For a discrete rotation with angle 7r/2 around the center, the transformations of global and 
relative coordinates are denoted by 



and that of fermion is 



N^R{N), r^R{r), 



(14) 



(15) 



V12 is a rotation matrix about a spinor index in the 5*0 (4) base, up to a phase factor given by 
a form 



V12 = ^F5(7i - 72) (1 + 7172) = e''' 



The following four operators Oi for i — 1,2, 3, 4, 



f -i \ 

i 

i 

\0 i / 



(16) 



Oi = 1, 02 = irn'2 = Ts, 03 = 71 + 72, O4 = r3(7i + 72), 



(17) 



are invariant under the rotation V\20iVi2- Our analyses in the following sections concentrate 
on the improved staggered fermion action by these four operators with U2N+a,fi = 1- 
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3 Analysis of Mass Matrices 



To split masses in degenerate tastes we introduce four rotationally invariant operators, which we 
denote as Aii = ^Oi'^ JT], for the original staggered fermion action Explicit expressions 
for Aii are given by 

r,r' 
r,r' 
r,r' 

M4{N) = J2^r{N){T,{^i+^2)}r,r''^r'{N). (18) 
r,r' 

The total mass matrix form which is invariant under the rotation by 7r/2 in two dimensions is 
given as 

Mr = mil + m2r3 + 7713(71 + 72) + m4r3(7i + 72) 

/ mi —im^ + 1714 —im^ — —im2 \ 



im^ + m4 mi —im2 —im^ + m4 

im3 — m4 2m2 mi im^ + m4 

V im2 im^ + m4 —im^ + m4 mi / 



(19) 



where mi,m2,m3,m4 are parameters of each operator in Eq. (fTS|). Mr has four eigenvalues 

mi — m2 — \/2m^ + \/2mi, mi — m2 + \/2m^ — A/2m4, 

mi + m2 — \/2m^ — -\/2^4, rni + m2 + \/2m^ + A/2m4. (20) 

The operator M.i cannot separate the 2- fold degeneracy in two-dimensional staggered Dirac 
fermions. The degeneracy can be solved by remained three operators M.21 and A 
4-component spinor should be separated into two 2-component spinors since a two-dimensional 
Dirac spinor is composed of a 2-component mode and we keep the rotational invariance even 
under a finite lattice constant*. Actually all possibilities of this separation are three cases and 
are listed in Table Q Let us analyze these mass matrices in order explicitly. 



*If one permits the rotational invariance only after taking the continuum limit, it is not necessary for 
degeneracy of a heavy mode and there are six more cases derived from Eq. (|20|l . 





parameter conditions 


rotationally invariant mass term 


mass eigenvalues 


case 1 


1712 = ms = 


Mri =mil + 771413(71 + 72) 


mi ± \^m4 


case 2 


1712 = 11^4 = 


Mr2 = 777i1 + 7773(71 + 72) 


mi ± \/2m3 


case 3 


777.3 = 1714 = 


Mrs = mil + ?772r3 


777.1 ± 7n2 



Table 1: Three cases for the mass sphtting into two spinors. 



Case 1 

The mass matrix 



M, 



Rl 



( mi 7774 —7774 \ 

7774 777 i 7T74 

—7774 777i 7774 

\ 7774 7774 777i / 



can be diagonalized as 



(21) 



M*"^ — pIMrxPi — diag{mi — V2m4, mi — V2m4, mi + V2m4, mi + V2m4), 



(22) 



by a unitary matrix 



2V2 



( l+i 1-7 1-Z 1+7 \ 

-y/2 -V2 V2 

\/27 -\/27 \/27 -\/27 

1-i 1+i 1+i 1-i 



(23) 



Note that the matrix Pi can diagonalize the rotation matrix V12 simultaneously. The unitary 



transformed spinor is given by 



2\/2 



/ 1 - 7 -v^ -v^i 1 + 7 \ 

1 + i -\f2 V2i 1-7 

1+7 ^2 -V2i I- 1 

V2 V27 1+7/1 ^4/ 

/ (1 - 7)^1 - ^^2 - ^27^3 + (1 + i)^4 \ 
(1 + 7)^1 - \/2^2 + ^i*3 + (1 - ^)^4 
(1 + 7)^1 + ^2*2 - V2«*3 + (1 - 0^4 

\ (1 - 7)*i + V2*2 + V^«*3 + (1 + 0*4 / 



/ *1 \ 

*2 
*3 



(TV). 



(TV) 



Case 2 

The mass matrix 



il2 



( mi —im^ 

ims mi 

ims 

V 77773 



-77773 




\ 

-im^ 
mi im^ 
-im^ mi ) 



(24) 



(25) 
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can be diagonalized as 



Mr2^ = P2MR2P2 = diag{mi - V2m3, mi - V2m3, mi + V2m3, mi + V2m3), 
by a unitary matrix 





( -^/2^ 


A/2i 




-72i\ 


1 


-l+i 


l+i 


l + i 


l-i 


2^/2 


-1 - i 


1 - i 


1 - i 


l + i 




\ 72 


72 


-72 


72 / 



with the transformed spinor 
Case 3 

The mass matrix 



272 



/ 72i^i-(l + i)^2- (1-0^3 + 72^4 \ 
-72i^i + (1 - i)^2 + (1 + i)^?, + 72^4 

72?:^i + (1 - i)^2 + (1 + ^)^3 - 72^4 

\^ 72«*i + (l + i)^2 + (1-^)^3 + 72*4 I 



/ mi —im2 \ 

mi — im2 

im2 mi 

V im2 TTli J 



can be diagonahzed as 

^Aofl ^ P^MrsP^ = diag{mi - m2, mi - m2, mi + m2, mi + m2), 
by a unitary matrix 



P3 

with the transformed spinor 

^Ms{N)^Pi^{N) = 



272 



272 



f -l+i l-i -V2i -V2i \ 

72i V2i l-i -l+i 

72 72 l+i -l-i 

\ l+i -l-i 72 72 y 



/ -(1 + i)*i - 72i*2 + 72*3 + (1 - «)*4 \ 
(1 + i)^i - 72i^2 + 72^3 - (1 - ?)*4 

72i^l + (1+ i)^2 + (1 - 0^3 + 72^4 

72i*i-(l + i)*2- (1-0*3 +72*4 J 



The feature of our formulation is to keep the discrete rotational invariance. After the 
splitting, we can find the character of a Dirac spinor under the rotation, 



V'(x) ^ Q^{R{x)), 







phase factor of V12 


case 


1 


Q \ 

[ e'-gt j 




case 


2 


( Q \ 

\ e^'^gt J 


gii? = e*^ = -1 


case 


3 


( Q' \ 

[ e^^/2(gt)2 ) 


= e-*-/4 = (1 - i)/v^ 



Table 2: The properties of Dirac spinors under the rotation. 



where Q = e 



(i7r/4)(T3 



gi7r/4 g 

g g-j7r/4 ) • Actually in cases 1 and 2 we can keep the property of 



a Dirac spinor on lattice, 



(34) 



By contrast, ^(A^) acts as a vector not as a spinor in case 3. The properties of 2-component 
spinors under the rotation are summarized in Table |21 



4 Pole Analysis and 2-point Functions 

Our adding terms do not commute with the staggered Dirac operator. As a result, our analysis 
in the previous section is insufficient to split masses. We must proceed in the pole analysis of 
the theory because a pole mass is physical. For the help, Dirac fields are Fourier transformed 
as 

<friN) = / -£-^,(p)e^P^, ^riN) = / -JL^^^pyvN^ (35) 



(27r)2 



and the action becomes 



Y: MN) (D;) (;v,^0 {T,,,) ^r,r')^r' {N') 
N,N',r,r',^i,T 



The staggered Dirac operator is explicitly written as 



^ {27^ sin + 27^(1 - cosp^)} 



^Ap). (36) 



Dst{p) = 5I{^7Msinp^ + i7^(l-cosp/,)} 



/ is2 + C2 isi + ci ^ 

is2 — C2 isi + ci 

isi — ci —is2 — C2 

\ isi - ci -is2 + C2 / 



(37) 



isi - ci -is2 + C2 
where Si = sin pi and q = 1 — cos pi for i = 1, 2, respectively. 

Our steps to find a pole mass are as follows: (i) set pi = and p2 = in (pure imaginary) of 
the inverse propagator in the momentum representation where our rotationally invariant 
operators (fTK|) are included; (ii) calculate four eigenvalues A of D^^; (iii) find values of k, in 
setting A = 0. Four values of k equal to pole masses. As mentioned in sections 2 and 3, 
we keep the rotational invariance in our action and generate two Dirac spinors with different 
masses. The light Dirac mass is denoted as mi and the heavy one is done as m/j. We note that 
parameters mi, m2 = —im,2, m'^ = —im^ and m^ are real to obtain real pole masses. 



Case 1 

The improved staggered Dirac operator is given by 

= {^7/^ s^^Pii + ^7m(1 - cosp^)} + mi + m4r3(7i + 72) 



/ mi is2 + C2 + m4 isi + ci — m^ \ 

is2 — C2 + m4 mi isi + Ci + m4 

isi — Ci — m4 mi —is2 — C2 + m,4 

isi — Ci + m4 —is2 + C2 + m,^^ mi / 



Setting pi = 0, we obtain eigenvalues of Eq. (pH|l as 



A = mi ± 



\ 



2m4 - 4sin^ ^ ± J -Ami ( sin^P2 + 4sin^ 



2 V -V ^ 2 
Taking p2 = ik aX \ = Q m. Eq. we have the equation for the pole mass 



16(1 - m4) sinh^ - - 8(mi + 2m4) sinh^ - + (mi - 2m4)^ = 0. 



(3^ 



(39) 



(40) 
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Hence we find pole masses for Im^l < 1 as 



±ln 



±ln 



2 + ml-VA 
2(1 - ml) 



2 + mf - v^' 
2(1 -ml) , 



2 + ml + y/A 
2(1 -ml) 



+ 



\ 



'2 + ml + y/A^ 
2(1 -ml) 



- 1 



±mi, 



= ±.mi 



(41) 



wliere A = m\ {%m\ + {m\ — 2m\)'^]. 

The mass is still splitting in the pole analysis for the improved staggered action. From 
Eq. 1)411) ■ it is found that we can take a limit \mh\ — > c>o for arbitrary m/ by performing e 
in an expression m\ 



^2-l-e(0<e<l). 



Case 2 

The improved staggered Dirac operator is given by 

= H {^7m sinp^ + z%(l - cosp^)} + mi + im3(7i + 72) 
= XI sinP/x + ^7m(1 + "^3 - cosp^)} + mi 



mi is2 + C2 isi + c'l 
— C2 mi 
isi — c'l mi 



v 



isi — c'l —is2 + C2 





isi + c'l 

-ZS2 - C'2 

mi 



where c[ = 1 + m'^ — cos Pi . Eigenvalues of Eq. (jl^ with pi = are 



2P2 



A = mi ± Y — 4(1 + m'3) sin' 

Setting p2 = in and A = 0, the pole mass is satisfied with 

o K m? + 2(m'n)'^ 

smh^ - = ^^4^- 

2 4(1 + m'3) 

Solutions of this equation under — 1 < mg are 

m ? -I- 

K = ±2 In 



ml + 2(m3)2 



+ 



m\ + 2(m^)2 



+ 1 



(42) 



(43) 



(44) 



(45) 



^ 4(1 + m^) ^ 4(l + m(5) 
The pole mass remains degenerate because the improved term m3(7i + 72) is absorbed into the 
kinetic term. 
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Case 3 

The improved staggered Dirac operator is given by 



m 

{ii7^ sinp^ + i7^(l - cosp^)} + mi + im'^T^ 
/ mi is2 + C2 isi + Ci 



«S2 - C2 



mi 



isi — ci — m2 
—m'2 isi — Ci 



m2 
mi 

-iS2 + C2 



isi + Ci 
-is2 - C2 
mi 



(46) 



Eigenvalues of Eq. (jiUI) with pi = are 

A : 



2P2 



nil i \/ — 4 sin — 



m 



/ \2 



± 4mf, sin^ — . 



Setting P2 = in and A = 0, the pole mass is satisfied with 



K 



I ^2 



sinh — 



m\ + {m'2) 



Ail ± m'o 



(47) 



(4J 



Solutions of this equation under — 1 < < 1 are 

±2 In 



K 



I ml + (m'a)^ 



4(1 + m2) 



+ 



/m^ + (m'g)^ 
4(1 + m2) 



+ 1 



±2 In 



mj + (m'g)^ /m^ + (m'g)^ 



+ 



1 



(49) 



4(1 - m'2) V 4(1 - m'2) 
This case allows pole masses to split although the rotational property of the eigenmode is not 
a spinor from the discussion of the previous section. 

In summary, we find that the improved staggered Dirac operator for cases 1 and 3 can 
split the degenerate mass through the analysis of the inverse propagator. However, that of 
case 2 cannot do because the additional effect is absorbed into the kinetic term of the original 
staggered Dirac operator. In addition, we should mention that there is no massless mode in 
cases 2 and 3 from Eqs. (|45p and (j49p . These nontrivial results are originated from the fact 
that the rotationally invariant mass terms do not commute with the staggered Dirac operator. 
Finally note that it is possible to take the light mass mi to zero by tuning mi and m^ only in 
case 1. Focusing on this particular situation, we discuss about massless and infinity modes in 
the next section. 
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5 Massless and Infinity Modes of Rotationally Invariant 
Action for Viable Even-Odd Separation 



In our two-dimensional lattice formulation, a matrix 

Ts = 7i727i72 = 



^10 0\ 

0-100 
0-10 

V 1 / 



(50) 



can define a chiral projection for a Dirac spinor because of {Dst, ^5} — 0. From the explicit 
matrix representation, the positive chiral mode is put on even sites and the negative mode 
corresponds to odd sites. 



( ^1 \ 




V *4 / 



1 - Tr, 



( \ 

V / 



(TV). (51) 



We note that the chiral projection is not discrete rotationally invariant since [^12, 7^ 0. 

This (even-odd) chiral property of staggered fermions holds in general D-dimensional cases 
First of all, we define chiral projection operators 



R 



(52) 



and a Dirac spinor ^ is projected out as 



The kinetic term of a staggered Dirac fermions lagrangian can be written as 



(53) 



(54) 



Among operators discussed in sections 3 and 4, A^3 and M.^ are chirally invariant but M.i and 
M.2 change chirahties of Dirac spinors. Therefore, M.i and M.2 can construct Dirac mass terms 
while and M.^ can do Majorana mass terms. Although we can define the formal discussion 
of the chirality just as normal chirality 75, this definition of the chirality depends on a special 
lattice frame because [V^,^, r2£i+i] 7^ where V^i, means the /xi^- rotation by 7r/2. 
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Without mention about the chirahty, we can define a massless mode and can throw up the 
mass of a heavy mode to infinity in case 1 of section 4. Solutions of Eq. ()40|) under the massless 
condition ml = 2m\ are determined as 



sinh 



smh 

2 



2 



0, 



2ml 



(55) 



(56) 



1 — ml' 

It must be noted that the eigenmode of the Dirac operator around a massless pole is not 
orthogonal to that around a heavy mass pole because their Dirac operators are different from 
each other. The massless modes for m^ > are explicitly written by 



f 1 + V2 \ 
-1-V2 
1 
1 



In m4 < case, they are expressed as 



f 1-V2\ 
1 - V2 
-1 
1 



(57) 



/ 1-V2 \ 
-1 + V2 
1 
1 



f l + V2\ 
1 + V2 
-1 
1 



In order to get the eigenmode for the heavy mass, we take pi 
D™^(pi = 0, P2 = irrih) is given by a form 



(5^ 



and p2 = imh for Eq. (fH^ . 



/ mi 1 + m4 — e" 

-1 + m4 + e"™'' 



V 



-m4 




mi 




— m4 


mi 
m4 — e" 





m4 

-1 +m4 + e™'' 
mi 



(59) 



m4 1 

Furthermore, from Eq. ()5fij) . we find that eigenmodes corresponding to the heavy mass are given 

by 

/ r ^ ^ J2ml{l + ml)] \ 

m4 I ^ J 



1 — m4 
m4 

— (1 — m4 



2m4 + m4 \/l + m^ 



2ml 



l + ml 



I — ml 



(60) 
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for rrih > and 



-(1 — m4 

1 — m4 



2ml 



m4 



2mj + m^Jl + ml 



V 



1 + m4 
m^ 



|m4(l + m^) + y'2m|(l + m|^ 



(61) 



for mh < 0. 

To decouple the heavy mode, we can throw the mass up to infinity. Actually from Eqs. (fS^ 
and we can realize massless and infinity modes as Table El simultaneously. From these 
results, it is found that infinity modes can be separately put on even or odd sites. 



7714 > 



m^ < 



massless modes 



/ 1 + V2 \ 
-1-V2 
1 
1 



/ 1 



V 



V2\ 
1 - V2 
-1 

1 / 



/ 1-V2 \ 
-1 + V2 
1 
1 



/ l + v^\ 
1 + V2 
-1 
1 



infinity modes 






V J 






V 1 / 



1 

1 ' 

V y V y 



Table 3: Eigenvectors of the improved Dirac operator in case 1 with mf = 2m\. Expressions 
of infinity modes are given by substituting m\ = 1 for Eqs. (j6Uj) and (j6ip . 



6 Summary and Discussion 

We have studied the mass splitting of two-dimensional staggered fermions based on the 5*0 (4) 
Clifford algebra. Introducing four rotationally invariant operators, we have analyzed three types 
of improved staggered Dirac operators and found one possibility (case 1) for taking a single 
mode in a two-dimensional free theory. The case keeps the splitting not only in the analysis 
of the mass matrix itself but also in the pole analysis including the kinetic term. According 
to the improvement with respect to the rotational invariance, the derived 2-component modes 
can be regarded as the ordinary spinor under the rotation by 7r/2. Furthermore, one can find a 
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massless mode in the case unexpectedly. The formal chiral projection which means even-site 
and odd-site separation of fermion modes is not consistent with the rotational invariance of a 
staggered Dirac action. Nevertheless, massless and infinity mode-representations in the case 
realize even-odd separation of the infinity mode. 

Our future tasks are analyses of interacting theories and the extension of our approach to 
four dimensions. In particular, it is crucial that the stability for the massless condition under 
quantum corrections by gauge interactions. Namely, in the case that the theory is not stable, 
it may be uninteresting that one needs a fine-tuning of the additional mass parameter as in 
Wilson fermions. For the infinity mode, it is very interesting if the even-odd separation is 
induced in staggered fermions even when we consider interaction effects. 
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